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Abstract. Recent developments in QCD phenomenology have spurred on several improved ap- 
proaches to Monte Carlo event generation, relative to the post-LEP state of the art. In this brief 
review, the emphasis is placed on approaches for 1) consistently merging fixed-order matrix ele- 
ment calculations with parton showers, 2) improving the parton shower algorithms themselves, and 
3) improving the description of the underlying event in hadron collisions. 
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INTRODUCTION 

The immediate horizon of accelerator-based high energy experiments is dominated by 
HERA, its legacy and final years of running at DESY, by the Tevatron, currently in 
its second run of operations at Fermilab, and by the Large Hadron Collider, under 
construction at CERN. A common denominator for all three machines is the study of 
^ ; high energy hadronic interactions at unprecedented levels of statistical precision. Thus, 

for a wide range of measurements, the limiting factors are ultimately systematic and 
theoretical in nature, rather than purely statistical. 
Q_i! Among the most important challenges is naturally that, while perturbative QCD de- 

scribes the interactions of quarks and gluons, experiments observe hadrons. In addition, 
many collider observables involve an interplay between widely separated energy scales, 
logarithms of which may appear and impact the validity of predictions even at the pertur- 
bative level. As a result, the field of QCD phenomenology is experiencing a rapid pace 
of development, a significant portion of which can be traced to either of two sources that 
I will focus on here. 

Firstly, as the Centre-of-Mass energy increases, the phase space for radiation also be- 
comes larger; high-,? final states are likely to be accompanied by high jet multiplicities. 
To accurately predict observables in such processes, lowest order scattering matrix el- 
ements are not sufficient. Rather, more sophisticated approaches are called for, which 
combine the rate of hard wide-angle jets predicted by fixed-order 2 — ► n matrix ele- 
ments with a resummation of multiple soft emissions in a consistent way, avoiding both 
double counting and "dead regions" over all of phase space. 

Secondly, hadron collisions involve new intrinsic challenges relative to ee (and to 
some extent also ep) scatterings, since both of the initial states are here composite and 
strongly interacting. However, on the theoretical side, the description of beam remnants 
and underlying events has gone through a long period of relative hibernation, essentially 
during the LEP era, with few new ideas emerging over the last ~ 20 years. Recently, 
however, interest has been rekindled, largely in response to increased interest from the 
Tevatron and LHC collaborations. 



Staying along the same lines, it should also here be emphasized that, when using 
LEP, HERA, and Tevatron results to make predictions for the LHC, an extrapolation is 
performed over orders of magnitude in Q 2 and x, at the same time as approximations 
that are "safe" at lower energies may be stretched into regions where large corrections 
are to be expected. As such, a non-trivial and many-faceted issue is how to treat the 
associated systematic and theoretical uncertainties. Though I will not touch directly on 
this topic below, some recent progress that partly addresses it is the emergence of parton 
distributions with intrinsic errors, reported on elsewhere in these proceedings JESS]- 



HARD & SOFT - ME/PS MATCHING 

The evaluation of tree-level transition amplitudes, involving less than, say, 5-6 partons 
in the final state, is a procedure which by now has been largely automated. Matrix 
Element Generators like CompHEP H, Grace [El 01, HELAS fl, MadGraph JHIl, 
O'Mega [ 10], and AMEGIC++ ifllll provide fast and reliable means of obtaining (more 
or less) tractable analytic expressions for a broad range of matrix elements, both in 
the Standard Model and beyond. Combining these with efficient numerical phase space 
integration algorithms such as BASES/SPRING [12] and others [.13], it is possible to 
further automate the phase space weighting, and hence to generate events corresponding 
to the chosen amplitude at matrix-element level. A number of more dedicated matrix 
element evaluation codes also exist, with processes hard-coded one by one, most notably 
AlpGen |QJ], but also MCFM H and others EES [IE El III SI. 

Common to all these approaches (most at tree level and a few at one loop) is that 
they represent fixed-order expansions in the electromagnetic, weak, and in particular 
strong coupling constants. As such, the virtue of these calculations is, briefly stated, 
that they include the entire helicity and interference structure of the amplitude (as well 
as virtual corrections to it, to the extent that loops are included), up to the given order 
calculated. Furthermore, asymptotic freedom implies that the stability of this expansion 
should improve with energy, due to the gradual vanishing of the strong coupling at 
large energies. Admittedly, the complexity rapidly grows with the number of particles 
involved, and so as already mentioned, these approaches are presently limited to fairly 
inclusive observables, where the number of resolved final state particles does not exceed 
a handful or so. 

On the other hand, from LEP we know that multiple soft gluon emissions are impor- 
tant in building up the full event structure. Mathematically, these corrections correspond 
to logarithmic enhancements of the amplitude, of the form 

«flog 2A '(ehard/esoft)+... (1) 

where <2soft/<2hard is a measure of the softness of the gluon(s) relative to the hard 
scale(s) in the problem. Moreover, when going to higher energies, the phase space for 
such emissions increases. Thus, while fixed-order calculations should be able to predict 
reliably the rate of a few well-separated jets (and other observables at the same level 
of inclusiveness), it is necessary go beyond the fixed-order approximation to obtain a 
picture of the full event structure. 



To improve the logarithmic accuracy, two dominant approaches exist: resummation 
calculations, and parton showers. Both are approximations to perturbation theory which 
work at infinite order in the coupling constant, and which are exact in certain limits. 

The former approach, resummation, allows to include not only the terms shown 
explicitly above (double logs ~ LLA), but also less singular logarithms in a systematic 
way. However, the formalism can still only be applied to relatively inclusive quantities, 
and a separate calculation must be performed for each observable, though interesting 
work has recently been carried out on automating calculations of this type ll22ll . 

In this talk, I will concentrate on the parton shower approach. While this description is 
formally correct only to leading logarithmic accuracy, it has the virtue that a fully exclu- 
sive description of the final state is obtained, which can be easily matched onto hadro- 
nisation descriptions, and from which in principle any observable can be constructed. 
Moreover, it is possible (and indeed necessary, as e.g. in the case of momentum conser- 
vation) to include at least a subset of higher-order effects, such as angular ordering of 
emissions, optimizing the scale choice in a s with respect to higher order kernels II23I1 . 
and choosing azimuthal angles in the branchings non-isotropically [24]. In practice, 
such refinements have been introduced in all of the standard shower Monte Carlos, in- 
cluding in particular the ARIADNE dipole shower EES, the HERWIGjSSjllc ] 
and HERWIG++ [31] showers, and both the PYTHIA virtuality-orderedJpJSH 
and transverse-momentum-ordered l35ll showers (the SHERPA generator [36] basically 
uses a variant of the PYTHIA virtuality-ordered algorithm). It would therefore be grossly 
misleading to equate leading log analytical calculations, where no such refinements are 
included, with leading log parton showers. 

The parton shower approximation starts from the observation that the collinear limit 
of QCD (and QED, for that matter) is universal. Thus, a process like e + e~ — > qq can 
be corrected to the process e + e~ — ► qqg using universal expressions for the q — > qg 
splitting probability. Due to the universality, the same expressions may then be applied 
again to describe the radiation of further gluons, as well as gluon splittings into quarks 
and so forth. Since the integrated probability at each step is nominally infinite, an 
ordering is introduced, whereby the emissions are generated sequentially according to 
some resolution criterion, like angle, virtuality, or transverse momentum. A lower cutoff 
on the resolution variable may then be naturally introduced, that regulates the infrared 
divergences, and at which scale a hadronisation description is supposed to take over. 

Thus, the virtues are that final states with an arbitrary number of partons may be 
built up, with a transition to hadronisation descriptions built in from the start. The down 
side is that the approximation is only exact in the collinear limit. For hard and/or wide- 
angle emissions, different parton showers can give widely different answers, reflecting 
the approximate nature of the approach in those regions. However, as mentioned above 
these are precisely the regions where the fixed-order calculations are at their best, hence 
it has been a long-standing wish to join consistently the state of the art of both worlds. 



TABLE 1. ME/PS Merging: List of processes, X, for which X+jet merging has been implemented 
in the HERWIG and PYTHIA Monte Carlos (e + e~ processes omitted). 





pp -> h° pp - 


-> V (= 7/Z/W) 


DIS 


top decay 


SM decays 


SUSY decays* 


HERWIG 


V 


V 


V 


V 






PYTHIA 1 


V 


V 


V 


V 


V 


V 



* corrections approximate for 2-body RPV modes, absent for RPV 3-body modes l37ll . 
' PYTHIA: applies to both the Q 2 and ordered shower algorithms in PYTHIA 6.3. 



ME/PS Merging 

The simplest (and oldest) approaches to join matrix elements and parton showers 
I will here refer to as matrix-element/parton-shower (ME/PS) "merging", to be con- 
trasted with "matching" below. Essentially, merging improves the parton shower off a 
hard system, call it X, by re-weighting the position of the hardest jet in phase space to 
reproduce the matrix element distribution for X+jet. An overview of hadron collider pro- 
cesses for which such corrections have been implemented in the HERWIG and PYTHIA 
models is given in Tab. [T] 

Technically, the way these corrections are implemented can be quite different, de- 
pending on the showering algorithm. In HERWIG, the showering algorithm has a "dead 
zone" in the hard wide-angle region, where no radiation at all is produced. In order to 
match to the matrix element which does produce jets there, two classes of events are 
effectively merged, e.g. X and X+jet, with the latter chosen inside the dead region of the 
former Il29ll . The detailed procedure is somewhat complicated and has only been worked 
out for a few cases, most recently for Higgs production H38h . 

In PYTHIA, the problem is rather the opposite. Too much radiation is generally 
produced in the hard wide-angle region, as compared to the matrix element answer. It is 
thereby straightforward to introduce a re-weighting, vetoing some of the extra emissions, 
to arrive back down at the matrix element rate 113 3l 13911 . Note that these corrections are 
applied for both the Q 2 - and -ordered shower algorithms in PYTHIA 6.3. 

So far, so good. However, for both the HERWIG and PYTHIA style merging, the 
procedure rapidly becomes more involved when attempting to generalize the methods to 
more complicated final states (see e.g. 114(1 141I1L Moreover, recent developments along 
related lines have resulted in a range of more generic approaches which are now being 
more actively pursued, as will be discussed below. 

ME/PS Matching at Leading Order 

The problem of consistently adding parton showers to a set of leading-order matrix 
elements for X, X + jet, X + 2 jets, etc, to obtain an inclusive sample of X production, 
matched to all available hard radiation matrix elements, has recently been studied in 
detail by a number of authors, in particular by Mangano (MLM) II42I1 . by Catani, Krauss, 
Kuhn, and Webber (CKKW) SQ], by Lonnblad |0], and most recently by Mrenna 
and Richardson 14611 . (See also talk by Frixione MM.) 



All these approaches essentially allow a consistent adding together of events gener- 
ated with different jet multiplicities at the matrix-element level (e.g. W, W+jet, W+2jets, 
...), by re-weighting them and showering them in way so that double counting and empty 
regions are avoided over all of phase space. 

The approach proposed by CKKW II43I1 is, briefly stated, to first select the jet multi- 
plicity, n, at the matrix element level according to a known probability, 

p ( n ) = ;7X77T77T77T — ' 5 = (2) 
00 + 0\ + 02 + 03 + . . . 

with <2 cu t a cutoff scale (in principle arbitrary) that regulates the infrared divergencies 
of the matrix elements. According to the chosen matrix element, a set of explicit four- 
momenta pi are then generated, to which a jet clustering algorithm is applied. A series 
of 'branchings' is thereby reconstructed, which can be interpreted as a parton shower 
history. The event as a whole is then re- weighted according to the Sudakov form factors 
(see below) and a s values associated with the reconstructed intermediate scales. A parton 
shower can then be applied as the final step, with emissions above the cut scale vetoed 
for all except the highest jet multiplicity matrix elements. 

Why this works is more technical: recall that the parton shower is formulated in terms 
of the no-emission probability between two scales, the Sudakov form factor, which 
in all simplicity is the (singular part of the) probability for an n-jet configuration to 
remain an n-jet one as a function of the resolution scale. By re-weighting the matrix 
elements with Sudakov factors for each leg, the leading divergencies which would lead to 
double-counting between e.g. 0\ and 02 are cancelled. Roughly speaking, the Sudakov 
re-weighting takes into account that for every 2-jet configuration you gain, you must 
lose one 1-jet one. It was shown by [14311 that this procedure makes the sum stable at least 
to next-to-leading logarithmic (NLL) accuracy. Note that this stabilisation is to some 
extent equivalent to the cancellation of real divergencies by virtual ones in a full NLO 
calculation, with the Sudakovs here playing the role of virtual corrections, that have the 
same structure (but opposite signs) as the tree-level divergencies. 

To further explain what the Sudakovs are doing, note that the parton shower is correct 
in the limit of strongly ordered emissions, i.e. in the limit that each successive scale 
is much smaller than the preceding one. In this case, the emission probability is large 
(i.e. the Sudakov = no-emission probability, is small), and the Sudakov re-weighting 
becomes quite important. At the other extreme are emissions which happen at similar 
scales; here, the Sudakovs are very close to one (the probability to go from one scale to 
another without emission becomes unity in the limit that the two scales coincide), and 
hence the matrix element results remain essentially unaltered here, as desired. 

This style of matching has since been implemented in the SHERPA event generator for 
several processes [36]. It was then noticed by Lonnblad W45I1 that a better matching could 
be obtained by replacing the analytical Sudakovs used in the original CKKW prescription 
by Sudakovs numerically generated by running the actual parton shower (so-called 
'pseudo-showers'). This is the approach implemented 114 811 in the ARIADNE generator 
ll26ll . Mrenna and Richardson subsequently made further refinements [46], applying the 
methodology also to hadronic collisions, using MadGraph and the HERWIG and PYTHIA 
generators, part of which work is stored in the PATRIOT event database [49]. 



TABLE 2. LO ME/PS Matching: List of processes, X, for which LO 
ME/PS matched generators/samples are available. 



e+e-^qq pp -> V (= y/Z/W) pp -> VV DIS 



ARIADNE y V* - V 

SHERPA y/ V y/ - 

PATRIOT - V - - 



* V = W only. 



Mangano's prescription ("MLM matching") 1E2I1 is similar but somewhat simpler in 
spirit than CKKW. In particular, it is based on clustering of events after showering 
and thus has a much simpler interface between the matrix element and parton shower 
generators. 

Tab.Elgives an overview of processes for which Leading Order matching is currently 
implemented/available . 



ME/PS Matching at NLO 

Above, the aim was essentially to describe real QCD radiation as accurately as pos- 
sible, over all regions of phase space. This was accomplished by consistently matching 
on parton shower descriptions of soft radiation to a set of tree-level matrix elements 
describing as many hard emissions as one cares to calculate, which in practice currently 
means up to 3-4 extra jets. 

However, since virtual corrections are not included, the normalization of the distri- 
butions is still only correct to leading order. Quite recently, the problem of matching 
parton showers to the full NLO theory, i.e. one-leg and one-loop corrections to the 
lowest order, has been addressed by several groups (see also talk by Frixione iP^ll '). 
Early approaches include the use of phase space slicing to separate resolvable and un- 
resolvable regions B5QL 15 ill , which, despite a number of initial successes, suffers from 
the drawback of not reproducing the perturbative expansion correctly. Important studies 
have also been carried out by the group though so far practical applications 

have been limited. 

Presently, the most mature NLO matching approach is the one put forth by Frixione, 
Nason, and Webber |55, 56j], which has been implemented in the program MC@NLO 
(essentially a superstructure built onto the HERWIG Monte Carlo). Another promising 
approach suggested by Kramer and Soper 1571 [58] has so far only been applied to 
e + e~ observables, though work is in progress to generalise it. A more complete list 
of processes is given in Tab.|3l 

Naturally, the big boon is that one automatically obtains cross sections normalised to 
NLO precision. A vice, as compared to the leading order matchings, has so far been that 
the NLO matrix elements only include the first (real) hard emission; subsequent emis- 
sions, even when hard, must still be generated by the parton shower, though recent work 
has been carried out on including also CKKW-style matching for higher jet multiplicities 
& 



TABLE 3. NLO ME/PS Matching: List of processes, X, for which unweighted NLO ME/PS matched 
event samples can be generated with the MC@NLO generator (weights = ±1). 



e+e -^qq pp— >h° pp — > V (= y/Z/W) pp — > VV pp — > QQ single top 
MC@NLO V \/ V V 1 " in progress 



* Note that, due to the simple structure of e + e , NLO matching can here effectively be obtained using the existing 
HERWIG and PYTHIA merging methods, and re-weighting by the loop factor (l+a s /n). 

* Q: heavy quark. 



A related issue is that, since MC@NLO is hard-wired to HERWIG, it is presently 
not possible to vary the parton shower model. Consider for instance the peak position 
of the Drell-Yan and h° p± spectra. This is where the bulk of the cross section sits, 
and this is, roughly speaking, the region that gets the most enhancement by the loop 
corrections. However, this region is also highly sensitive to multiple soft emissions, 
which are resummed by the parton shower. A much awaited future development is thus 
the creation of tools that allow a more generic interfacing with different shower models. 

In a similar vein, note that, while the overall normalisation is improved by going from 
LO to NLO matching, the same is not necessarily true for the shapes of distributions, as 
follows. Consider again the Drell-Yan p± spectrum in hadron collisions. As mentioned 
above, it is dominated by Z + multiple soft emissions in the peak, while Z+jet dominates 
in the tail. It would now be possible experimentally to define a semi-inclusive Z+jet 
fraction as a function of p±z- At NLO, the prediction for Z+jet (and hence for this 
subsample) is effectively only correct to leading order, since virtual corrections to Z+jet 
are not included. Since the Z+jet fraction tends towards unity at large p±z values, the 
overall shape should probably not be regarded as being more precisely determined here 
than in the leading order matching schemes discussed above. The extreme example 
would be an observable sensitive to multiple hard emissions, where a leading order 
matching including several hard jet radiations would clearly be superior to the present 
NLO matching schemes, where the second jet has to be radiated by the parton shower. 



NEW PARTON SHOWER ALGORITHMS 



Both HERWIG++ [60] and PYTHIA 6.3 (611] contain new parton shower models. In the 



HERWIG++ case, refinements have been made IBlll on the already existing HERWIG 
model, while a completely new shower model IB5I1 has been implemented in PYTHIA 
6.3. The recently completed program APACIC++ 16211 also contains a parton shower 
model, which essentially is an adaption of the old PYTHIA shower model with the 
specific implementation of CKKW style matching in mind. 

The HERWIG shower is based on a strictly angular-ordered sequence of emissions 
lEll . This correctly accounts for coherence effects in the emission of soft gluons, but 
has the disadvantage that it leaves a 'dead zone' in the hard 3— jet region, which has to be 
filled in separately, as discussed above. The shower evolution is stopped once a fixed low 
scale is reached, at which time a transition is made to a non-perturbative hadronisation 
description, the HERWIG one being based on the cluster model [63]. 



P±evo\ — 

z(\-z){ml-ml) 

FIGURE 1. Final-state branching in the /?x-° r dered shower. The evolution is performed on each side 
of the dipole separately, in the variable Pj_~,i, here expressed in terms of the parton virtuality, m*, its rest 
mass, niQ (for massive partons), and the energy sharing fraction z appearing in the splitting kernels P(z). 
Kinematics are constructed in the dipole CM frame, conserving energy and momentum inside the system. 



The HERWIG++ algorithm JUl] starts from the same basic principles and thus inherits 
most features from the old shower, but the definition of the evolution variable has been 
slightly changed, to allow a more correct treatment of heavy quark radiation as well 
as a more consistent behaviour in the soft-gluon limit. Finally, the treatment of cascade 
decays of resonances interspersed with showers should be improved by a more consistent 
implementation of multi-scale showers. Recent studies exploring this algorithm can be 
found in ||3ll 165|] . At present, HERWIG++ is still mainly a tool for e+e collisions, 
though work is underway to extend it to hadron collisions. Another item on the agenda 
is the inclusion of a CKKW- style ME/PS matching scheme. 

More radical changes have occurred in PYTHIA 6.3, where a new transverse- 
momentum-ordered shower based on a recently proposed hybrid between the dipole 
and parton shower formalisms [35] has been implemented (in addition to the old 
virtuality-ordered shower, which is carried over from previous versions). The choice of 
p± as evolution variable was made for several reasons. Firstly, it has the dual property 
of simultaneously being a good measure of hardness while still leading to a natural 
angular ordering of emissions 116611 . Secondly, it is Lorentz invariant under longitudinal 
boosts, which is not the case e.g. for the HERWIG angular-ordered evolution variable. 
Thirdly, while both PYTHIA and HERWIG can be tuned to give good descriptions of 
the LEP data, the ARIADNE p^-ordered dipole shower tends to do even better. Note, 
however, that the p± measure proposed in l35ll is different from the one used in the 
ARIADNE evolution. Fourthly, the underlying-event model in PYTHIA is based on a 
/^-ordered sequence of multiple perturbative interactions, hence a shower also ordered 
in p± allows for a more unified treatment of the perturbative activity as a whole, as will 
be discussed below. 

The dipole/parton shower hybrid approach implies an evolution of individual parton 
lines, but with recoils occurring inside dipoles, as illustrated in Fig.[T] For the final-state 
shower, these dipoles are normally defined by the respective colour partners, with some 
exceptions involving decays of heavy resonances. For the initial-state shower, where a 
backwards evolution is performed from the hard scattering down, only a single dipole is 
relevant, spanned between the two incoming partons. 



Tevatron: ttbar + 1 jet ctfosi , nn K-fantnrs 
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FIGURE 2. The Tevatron tt+jet rate as a function of jet p±. Results are shown for MadGraph, and for 
hard ('power') and soft ('wimpy') variants of both the Q 2 - and /^-ordered PYTHIA showers, see text. 



Studies carried out so far indicate that the new shower leads to an improved de- 
scription e.g. of event shapes at LEP [67] and of Drell-Yan production at the Tevatron 
SH. Fig. ei gives a preliminary comparison of the tt+jet rate at the Tevatron, as a 
function of jet p\_ [69]. Results are shown for MadGraph (thick black line), for two vari- 
ants of the old <2 2 -ordered shower (green lines, solid and dotted), and for two variants 
of the new p^-ordered shower (blue lines, solid and dotted). The 'power' and 'wimpy' 
versions of each shower represent a reasonable range of variation of the maximum scale 
for shower emissions, with the former populating all available phase space, and the lat- 
ter bounded by the scale of the hard interaction. Though it is still too early to draw 
more general conclusions, it is encouraging that both the shape and normalisation of 
the distribution appear to be improved, especially in the 'power shower' variant of the 
p^-ordered algorithm. 



THE UNDERLYING EVENT 

The underlying event may be defined somewhat loosely as the activity in a (single) 
hadron-hadron collision which does not originate directly from the hard scattering 
that triggered the event. Currently, this component is not well understood from first 
principles. What is known, on the other hand, is that it produces omnipresent systematic 
as well as random fluctuations in activity, which impact isolation criteria, jet energy 
scales, etc., which can be of significant magnitude for experimental analyses. 

During the last few years it has become increasingly clear that the underlying event 
not only contains soft activity but also has a semi-hard, 'lumpy' component Il70ll . To 
explain this, the concept of multiple perturbative interactions [71] appears increasingly 
attractive. Several recent implementations are built on ideas incorporating multiple 
interactions in one form or another, including the underlying-event framework in the 



HERWIG add-on JIMMY E], the new interleaved model Jpjzl in PYTHIA 6.3, and 
the underlying-event model being developed for SHERPA Il74ll . this latter being quite 
similar to the old PYTHIA one 0. 

The most recent development has been the interleaved model of In this pic- 

ture, the evolution of the initial state radiation and the generation of multiple scatterings 
are no longer independent. Instead, a successive 'fine-graining' of all the perturbative 
activity is performed, with radiations interspersed (or 'interleaved') with interactions. 
This allows correlations to be introduced between all the perturbative activity at succes- 
sively smaller scales. For instance, the third interaction will know about the presence of 
a gluon having been radiated off a valence quark in the first, etcetera. This introduces 
for the first time non-trivial correlations in flavour and x, while all known sum rules are 
still respected (e.g. momentum conservation and quark counting rules). In addition, the 
model allows new possibilities for impact parameter dependence, and contains a more 
refined treatment of beam remnants, based on an extension of the Lund string model to 
'baryonic' colour topologies 1 3711 . 

With the development of more sophisticated physics models, the hope is that it will be 
possible to ask a range of more meaningful physics questions, among which especially 
the topic of colour correlations and colour reconnections (in view of a possible difference 
between the vacua left by e + e~ and hadronic initial states) is presently the most actively 
investigated. Finally, the energy dependence of the underlying activity is currently very 
poorly known; further studies aiming to pin down the scaling behaviour, for instance by 
including the available RHIC pp data, would be of great interest. 



CONCLUSION 

The quest for the cause of electroweak symmetry breaking is (hopefully) nearing its end. 
Should the Tevatron not discover it within a few years, the experimental programme at 
the LHC promises a comprehensive exploration of the TeV scale, with the emphasis on 
precisely this question, as well as on more general searches for signals of physics beyond 
the Standard Model. In addition, the high energies and hadronic environment of both the 
Tevatron and in particular the LHC will challenge our understanding of QCD, especially 
our ability to solve it for large numbers of partons, both real and virtual, and our control 
of phenomena near the borders of its perturbative domain. 

It is encouraging, then, to see a truly impressive effort mounted in recent years to 
address several of the most important issues, among which I have here discussed more 
precise and realistic simulation of events at higher perturbative orders, improved parton 
shower models, and some progress towards a better understanding of the underlying 
event. 

On a general note, the emergence of many sophisticated but specialised tools entails an 
increased need for efficient cross-communication between the programs. For instance, 
while the traditional Monte Carlo generators have been largely self-contained in the past, 
the time-consuming process of hard-coding matrix elements for each process by hand 
can now be left to automatic programs optimised specifically for this task. One should 
therefore expect an increased reliance on external interfaces, such as provided by the 



Les Houches Accords in the future. 

Finally, though much of the effort reported here is still centred around a relatively few 
groups, for which simple manpower restrictions often represent a non-trivial problem, 
the many new and exciting developments do seem to have led to an increased com- 
munication, bringing people together from different fields. Hopefully, with continued 
nourishment, this is a trend that will continue and grow in the future. 
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